Stochastic resonance with matched filtering 
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Along witli the development of interferometric gravitational wave detector, we enter into an epoch 
of gravitational wave astronomy, which will open a brand new window for astrophysics to observe 
our universe. Almost all of the data analysis methods in gravitational wave detection are based on 
matched filtering. Gravitational wave detection is a typical example of weak signal detection, and 
this weak signal is buried in strong instrument noise. So it seems attractable if we can take advantage 
of stochastic resonance. But unfortunately, almost all of the stochastic resonance theory is based 
on Fourier transformation and has no relation to matched filtering. In this paper we try to relate 
stochastic resonance to matched filtering. Our results show that stochastic resonance can indeed be 
combined with matched filtering for both periodic and non-periodic input signal. This encouraging 
result will be the first step to apply stochastic resonance to matched filtering in gravitational wave 
detection. In addition, based on matched filtering, we firstly proposed a novel measurement method 
for stochastic resonance which is valid for both periodic and non-periodic driven signal. 

PACS numbers: 04.80.Nn, 05.40.-a 



I. INTRODUCTION 



Since the concept of stochastic resonance (SR) was 
originally put forward in the seminal papers by Benzi 
and others ^ , it has continuously attracted considerable 
attention SR was subsequently found to occur in 
many kinds of nonlinear systems , even linear systems 
with multiplicative noise |4'| and diverse fields. Not only 
classical system but also quantum system [5] show SR 
behavior. Generally, SR emerged as a paradigm whose 
universal character was shown to be intimately related to 
the pervasiveness @ of the fluctuation-dissipation theo- 
rem Q . The possible use of SR in connection with weak 
signal detection experiments [sl-fisj. with special refer- 
ence to gravitational wave detection, was suggested since 
the infancy of SR . The application of stochastic res- 
onance to gravitational wave detection includes two as- 
pects, one is data analysis Q, and the other is related to 
detector itself fld\. As to data analysis method, the re- 
search about stochastic resonance almost all use Fourier 
transformation method to deal with the data with noise. 
On the other hand, the data analysis method in gravita- 
tional wave detection community is mainly the so called 
matched filtering method, i.e. maximizing the likelihood 
ratio over the parameters of expected signal templates 
So it is interesting to investigate whether we 
can combine the advanced data analysis method used 
in gravitational wave detection, matched filtering, with 
stochastic resonance. In this work, we will study whether 
the stochastic resonance phenomena can be observed also 
with matched filtering method to do data analysis. Our 
work shows that stochastic resonance does emerge with 
matched filtering data analysis method for periodic sig- 



nal, corresponding to the periodic gravitational waves; 
chirp signal, corresponding to the gravitational wave ra- 
diated from merging black holes or neutron stars and 
even random pulse signals. 

Till now we have three kinds of methods to character- 
ize stochastic resonance. The first one is power spectrum 
include output signal power itself and the ratio of output 
signal power to noise background. This kind of measure- 
ment method is most convenient and popular when deal- 
ing with stochastic resonance. But the disadvantage of 
this method is that it can only deal with periodic stochas- 
tic resonance. The second method is residence time dis- 
tribution [20j . The third one is information theoretic 
measures which includes the Shannon transinformation 
rate 21 1; the Fisher information 1221 : Kullback, Shan- 
non and other kinds of entropies [23ll: and information- 
theoretic distance measures [2J,[2^. In this aspect, our 
work shows that matched filtering method can be looked 
as a novel measurement method to characterize stochas- 
tic resonance. The method of matched filtering is very 
robust that it can characterize both periodic and non- 
periodic stochastic resonance. 

This paper is organized as following. In section II, we 
will give a brief review of matched filtering method used 
in gravitational wave data analysis. Then we present out 
our model equations used in this paper in section III. In 
section IV, we give out our result on stochastic resonance 
with matched filtering method. At last we discuss our 
result and its implications in section V. 



II. BRIEF REVIEW OF MATCHED FILTERING 
TECHNIQUES USED IN GRAVITATIONAL 
WAVE DETECTION 
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First we flesh out the brief description of matched fil- 
tering techniques used in gravitational wave detection 
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Assume our data stream takes form like 



h{t) = s{t)+n{t), 



(1) 



where s{t) is signal, while n{t) is Gaussian white noise. 
The likelihood ratio is defined as [In UM 
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Here u is the template determined by parameter in some 
space used to filter the data stream h{t). The maximal 
value is taken with respect to the time delay tc of the 
template. When the signal is periodic or almost periodic 
the maximal value can be equivalently taken respect to 
the initial phase of template. And the inner product is 
given by 
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ft is the noise-weighted cross-correlation between a and 
b. Here a(/) stands for the Fourier transform of a(t) 



a{t)dt, 



(4) 



and * stands for complex conjugate. Sn{f) is the one- 
sided power spectral noise density of the noise stream 
n{t) 



5„(/)EE2E[|n(/)|^ 



(5) 



where E[ ] denotes the expectation value over an ensemble 
of realizations of the noise. Since our data are all real, 
we can use one-sided power spectral instead of two-sided 
one. Above description is in frequency domain. In time 
domain, the inner product ([3]) takes form [13] 



{a,b) = 



a(t)'w{t — T)b{T)dtdT, 



(6) 
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where w{t) is Wiener's optimal filter [2a], which can be 
approximated as 
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Then the task of data analysis in gravitational wave de- 
tection is to maximize the likelihood ratio A4 over the 
parameters space of expected signal templates and then 
fix a detection threshold on this statistics. The detection 
threshold is simply established by the desired level of 
false alarm probability. The maximum likelihood is com- 
pared with the aforementioned threshold and a detection 
is announced whenever the threshold is crossed. 

In most pragmatic cases, the assumption of Gaussian- 
ness and whiteness to the noise is not valid, especially 



when noise passes through a nonlinear system, such as 
the system considered in this work. So the Wiener's fil- 
ter ([7]) is not valid any more. Therefore, we adopt the 
matched filtering method and the idea of Horner effi- 
ciency in optics |27J to define matched ratio as 
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where h{t) is the data stream and u{t) is the template 
as above, tc is the time delay for the template. A is 
the time delay range where the delay time is taken from. 
Theoretically A = (— oo,-|-oo), but it is impossible to 
take this infinity range practically. In practice, if the 
system has some characteristic time scale (such as the 
period of periodic system), this A can be taken as this 
time scale directly. Otherwise, there is an ambiguity to 
set this range. In the following sections we will come 
back to this problem and propose one trick to solve it. 
Different to equation (1), now n{t) may be any kind of 
noise, white or colorful, Gaussian or non-Gaussian. So 
our quantity will be valid in more general case than 
equation (2). Here the template is the signal profile we 
expected, and the different expected signal profiles up 
to a scale are looked as one template because we have 
already taken normalization in equation ([S]) . Since we do 
not know when the expected signal will begin in the data 
stream, we need take the maximal value respect to the 
delay time tc. 



III. MODEL EQUATION 

The model equation used in this paper is the famous 
overdamped bistable system [28| 

X ^aX-bX^ + Sp{t) + T(t), (9) 



with 



{r{t))^o,{r{t)r{t')}^2D6{t~t'). (lo) 



a and b are two parameters. We set a = 6 = 1 in this pa- 
per. This system has been widely investigated in stochas- 
tic resonance literature (28l - l30| . But all the studies used 
Fourier transformation method to do the data analysis. 
While in this paper we will adopt popular method used in 
gravitational wave detection — matched filtering method 
instead. Sp{t) is the signal which is the input signal of our 
system. If we take our model equation as the toy model 
equation of gravitational wave detector, Sp{t) is nothing 
but the gravitational wave signal. If we take our model 
equation as the afterward analysis tools for gravitational 
wave data analysis, Sp(t) is the raw data stream of the 
gravitational wave detector's output, which includes sig- 
nal itself and noise. In the following we take the for- 
mer viewpoint. As to numerical method, the stochastic 
Runge-Kutta algorithm |3l| is used in this work to sim- 
ulate equation The time step is set as At = 0.1. 
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We take 200 systems to realize the random process 
and take the ensemble average to analyze the random 
data stream produced by numerical method. Then our 
measurement quantity of stochastic resonance becomes 
the ensemble average of equation ([8]), i.e. i?[A^]. More 
explicitly, the measurement quantity we introduced in 
this work to characterize stochastic resonance is 



N 



max ■ 
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Here hi{t) is the data stream from i'th system. N is 
the number of the systems in the ensemble. In prac- 
tice, h and u are two arrays of numbers. Let us as- 
sume the length of the array is m. Then the convolution 
/-oo hilt — t)u(t — tc)dT can be done with standard algo- 
rithm. And we adopted rectangular method to approxi- 
mate the other integral. We need pay more attention to 
the time delay problem since our data array is finite. In 
practice, we take the elements from number to to of 
original discretized data stream u for delayed template, 
while the data stream h is correspondingly taken as the 
elements from number 1 to m — -f 1 of original array. 
Here A< is the sample time between the nearest two ele- 
ments of the array which is the time step of the numerical 
simulation in our work. 



expected profile of output signal we take Spit) directly 
as our template. In numerical simulation we evolve 6000 
steps which equivalents to 600 units of time which cor- 
responds to about 20 periods of input signal. Since the 
template here is periodic function, we can set the range 
A for the delay time tc as (0, The numerical result 
are presented in FigHJ In this figure we plot the i?[A^] 
(solid line) together with output signal power (dashed 
line) measured with traditional Fourier transformation 
method. The two results are consistent to each other 
very well. 
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IV. STOCHASTIC RESONANCE WITH 
MATCHED FILTERING METHOD 

In this section we present our main result on stochastic 
resonance with matched filtering method. We find that 
stochastic resonance shows up perfectly with character- 
istic quantity (|lip no mater what kind of input signal we 
use. Our input signal includes periodic one, which corre- 
sponds to the periodic gravitational waves; chirp signal, 
corresponding to the gravitational wave radiated from 
merging black holes or neutron stars and random pulse 
ones. All of these stochastic resonance phenomena can 
be described well by the quantity (|lip we first introduced 
in this paper. 



A. periodic input 

Firstly, we set the input signal in equation ^ as sinu- 
soidal function of time 

Sp{t) = Acos{ujt), (12) 

with uj = 0.2, A = 0.05 closely following [28]. This signal 
can model periodic gravitational wave very well. Peri- 
odic gravitational waves are emitted by various astro- 
physical sources. They carry important information on 
their sources (e.g., spinning neutron stars, accreting neu- 
tron stars) and also on fundamental physics, since their 
nature can test the model of general relativity. As to the 



FIG. 1: (color online) Stochastic resonance behavior of over- 
damped bistable system ([9]) for periodic signal with matched 
filtering measurement method. Ensemble averaged matched 
ratio(i7[A^] Ea. (|ll[) . solid line) and the output signal power 
(dashed line) are plotted. Both show perfect stochastic reso- 
nance behavior. In addition this two results are consistent to 
each other very well. The line marked with triangle is 
without time delay. We can see that i5[A1] with and without 
time delay show the same qualitative behavior. This is be- 
cause that our template actually has no time delay respect to 
the output data stream. 



B. chirp signal input 

Secondly, we set the input signal as the gravitational 
wave signal radiated from two merging black holes or 
neutron stars which reads 

Sp{t) = ^ cos [u;(tcoai - t)f"' . (13) 

A and uj are two parameters determined by the masses 
of the binary system's sub-bodies and the distance be- 
tween the source and the detector, tcoai is the time of 
coalescence. This wave form is theoretical prediction 
of post-Newtonian analysis for binary systems which in- 
cludes black hole-black hole binary, black hole-neutron 
star binary and neutron star-neutron star binary. In this 
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work we set A = 0.2, tcoai = 601 and uj = 6. This sig- 
nal is a typical wideband signal. Our parameters setting 
makes the frequency band locate around 0.2. The same 
to above subsection, the simulation time is 600 time units 
which is very near the coalescence time 601. 

Since this signal is not periodic, we have no time scale 
to guide us to choose the range of time delay A in equa- 
tion (jlip. Here we propose to use try error method to 
find out the true stochastic resonance behavior. Try er- 
ror method we mean trying different time delay range 
from small to large till the result converges. Here con- 
verge means resulting in the almost same optimal noise 
intensity. Theoretically the time delay is longer the con- 
vergence is better. But the data stream has finite length 
in practice. These results are plotted in Figl2] Since our 
input signal and template has no time delay at all, we ex- 
pect the line without time delay should be the converged 
result. But we see a small but nonvanishing difference 
between different time delays. This is because our data 
stream is some short, only 600 steps, while different time 
delay results in different effective length of data array. 
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FIG. 2: (color online) Stochastic resonance behavior of over- 
damped bistable system ([9]) for chirp signal with matched 
filtering measurement method. The ^^[Al] without time de- 
lay (marked with square) and with different time delay (10 
marked with disk, 20 marked with up triangle and 30 marked 
with down triangle) are plotted respectively. All of these lines 
show well convergence behavior which implies that the one 
without time delay has already described the true stochastic 
resonance behavior already. This is consistent with expected 
result, because our input signal and template has no time 
delay at all with equation ([9]). 



C. random pulse signal 

At last but not at least, we use random pulse signal as 
the input. We follow [2^ closely to construct the random 



pulse. Each pulse takes a value 

Sp{t) = A or Sp{t) = -A 



(14) 



and persists for time ^. Each A or —A are randomly 
taken with probability 1/2 respectively. In this work we 
set A = 0.05 and T = We set simulation time 600 
units still. 

Similar to above chirp case, we use try error method to 
determine the stochastic resonance again. We concluded 
all these results in FigEl 
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FIG. 3: (color online) Stochastic resonance behavior of over- 
damped bistable system ^ for random pulse input with 
matched filtering measurement method. This result is very 
similar with above ones but shows difi'erent optimal noise in- 
tensity. This is reasonable because of the difi'erent driven 
signal. This figure shows that matched filtering method can 
also describe non-periodic stochastic resonance very well. 

All three figures show that the one without time de- 
lay has already described the true stochastic resonance 
behavior already. This is consistent with expected re- 
sult, because our input signal and template has no time 
delay at all with equation ([9]). But generally, the time 
delay is important for the stochastic resonance measure- 
ment quantity (|lip when the template does has time de- 
lay compared with the data stream. 



V. CONCLUSION AND DISCUSSION 

Since 2003 the ground-based interferometric gravita- 
tional wave detector Laser Interferometer Gravitational 
Wave Observatory (LIGO) has taken several real exper- 
imental data. And other high sensitive detector such as 
VIRGO, GEO and TAMA have also achieved their ex- 
pected sensitivities. At the same time, the space-based 
interferometric gravitational wave detector LISA is also 
under way. Along with the development of these high 
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sensitive detector, we are entering into an epoch of gravi- 
tational wave astronomy. Gravitational wave detection is 
a typical example of weak signal detection where stochas- 
tic resonance may be applicable. But almost all of the 
data analysis method in gravitational wave detection are 
based on matched filtering. In contrast, almost all of the 
stochastic resonance theory is related to Fourier transfor- 
mation which does not relate to matched filtering. Here 
our work do the first step trying to relate matched filter- 
ing with stochastic resonance. 

The encouraging result of our work is that we do find 
that the stochastic phenomena emerges with matched 
filtering data analysis method. And more interestingly, 
matched filtering method can deal with not only periodic 
but also non-periodic stochastic resonance. Specifically, 
we show that matched filtering can deal with stochas- 
tic resonance which is driven by sinusoidal signal, which 
corresponds to periodic gravitational wave; chirp signal, 
which corresponds to gravitational wave from merging 
binary system; and random pulse signal. The interest- 
ing problem in the following is whether and how to ap- 



ply our findings and stochastic resonance theory to do 
gravitational wave data analysis and gravitational wave 
detector construction. 

On the other hand the measurement quantity (|1I|) , first 
proposed in this work, is a novel characteristic measure- 
ment for stochastic resonance. Our quantity can not only 
deal with periodic but also non-periodic stochastic res- 
onance, which improves the conventional Fourier trans- 
form method. We hope our new stochastic resonance 
measurement can find its usage in the future research. 
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